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ABSTRACT 


This  technical  memorandum  presents  the  results  of  analytical,  numerical,  and  experimental 
studies  predicting  the  velocity  and  position  of  a  high-speed  projectile,  Velocity  and  position 
predictions  were  made  using  computer  modeling  and  simulation  to  find  a  reliable  and  cost- 
effective  velocity  sensor.  The  two  most  popular  magnetic  field  sensing  devices,  the  copper  wire 
sense  coil  and  the  Hall  effect  sensor,  were  examined.  An  axisymmetric,  2-D  finite  element 
analysis  model  was  developed  (1)  to  calculate  the  field  distribution  due  to  a  permanent  magnet 
embedded  in  a  projectile  and  (2)  to  verify  the  analytical  model.  Experiments  were  conducted  and 
confirmed  the  accuracy  of  both  the  analytical  and  finite  element  models. 
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1.  INTRODUCTION 


The  main  focus  of  this  particular  research  project  was  (1)  to  predict  the  velocity  and 
position  of  a  high-speed  projectile  that  contains  an  embedded  rare-earth  permanent  magnet  (PM) 
using  computer  modeling  and  simulation  and  (2)  to  find  a  feasible,  reliable,  compact,  and  cost- 
effective  electromagnetic  motion-detection  sensor. 

This  technical  memorandum  (TM)  concentrates  on  the  first  aspect  of  the  research: 
developing  an  analytical  model  that  can  predict  velocity  and  position  of  a  high-speed  projectile. 
The  application  of  and  results  obtained  from  two  popular  magnetic  field  sensing  devices,  the 
copper  wire  sense  and  the  Hall  effect  sensor,  are  examined.  Because  these  sensors  must 
accommodate  uncertainty  in  the  actual  trajectory  of  the  projectile,  the  sensing  elements  must  be 
located  a  significant  distance  from  the  ideal  projectile  path.  This  required  distance  reduces  the 
magnetic  field  strength  observed  at  the  sensor,  thus  significantly  complicating  the  measurement. 

An  analytical  model  was  developed  to  predict  the  characteristics  of  the  voltage  waveform 
induced  in  a  sense  coil.  That  model  also  helped  define  the  requirements  of  the  data  acquisition 
and  peripheral  instrumentation  control  electronics. 

A  finite  element  analysis  (FEA)  model  was  also  developed  to  calculate  the  field 
distribution  due  to  the  PM  embedded  in  the  projectile  and  to  verify  the  analytical  model. 
Experiments  were  conducted  using  a  low-speed  (40  m/s)  projectile.  These  measurements 
confirmed  the  accuracy  of  the  analytical  and  FEA  models. 


2.  COMPARISON  OF  SENSING  DEVICES 


The  magnetic  field  sensor  appropriate  for  a  particular  application  depends  on  the 
performance,  noise  immunity,  reliability,  thermal  management,  economics,  and  packaging  of  the 
device.  Sense  coils  (figure  1)  and  Hall  effect  sensors  (figure  2)  are  the  two  most  common 
methods  of  sensing  the  magnetic  field  of  a  moving  PM  and  measuring  the  velocity  and  position  of 
that  PM.  In  most  motion  control  applications,  the  sensor  is  part  of  a  closed  loop  control  system. 
In  such  an  application,  the  sensor  is  fixed  and  external  to  the  system:  it  simply  measures  the 
magnetic  field  due  to  a  moving  PM  from  a  fixed  distance. 

The  voltage  induced  in  a  simple  sense  coil  is  a  function  of  the  magnetic  flux  present  at 
right  angles  to  the  conducting  coil  and  the  frequency  of  the  time-varying  magnetic  field.  This 
relationship  is  expressed  as 


Figure  1.  Sense  Coil  with  a  Moving  PM 


Figure  2.  Electron  Movement  in  a  Hall  Effect  Device 
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where  Fis  the  generated  voltage  at  the  coil,  N  is  the  total  number  of  turns  in  the  coil,  —  is  the 

dt 

total  derivative  operator,  and  <j>  is  the  magnetic  flux  exposed  to  the  coil. 

The  Hall  effect,  discovered  by  E.  F.  Hall  in  1879,  is  the  basis  for  all  Hall  effect  devices  [1], 
(A  conceptual  model  for  the  basic  principle  of  Hall  effect  sensor  is  shown  in  figure  2.)  When  a 
magnetic  field  B  is  applied  perpendicularly  to  a  current  Iref  flowing  in  a  semiconductor  bar,  an 
electric  field  is  produced  perpendicular  to  both  B  and  Iref  (proportional  to  BIref)  and  an  open- 
circuit  voltage  VH  is  observed.  The  ratio  VhJV,  where  F  is  the  applied  voltage,  is  proportional  to 
HHB1 ,  where  fiH  is  the  Hall  mobility.  When  this  physical  effect  is  combined  with  modem 
integrated  circuit  (IC)  technology,  many  useful  magnetic  sensing  devices  are  possible.  The  Hall 
element,  when  properly  biased,  produces  an  output  voltage  that  is  proportional  to  a  magnetic 
field.  To  make  a  linear  magnetic  sensor,  this  small  voltage  is  processed  through  a  high-quality 
amplifier,  which  produces  an  analog  output  that  is  proportional  to  the  applied  flux  density. 

The  Hall  voltage  is  defined  by 

where  VH  is  the  Hall  voltage  in  volts,  Rh  is  the  Hall  coefficient  (device  material  dependent),  d  is 
the  device  thickness,  Imf  is  the  current  across  the  device  in  amperes,  and  B  is  the  magnetic  flux 
density  incident  on  the  device  in  gauss. 

For  the  same  current,  electrons  must  travel  faster  in  a  poor  conductor  than  they  must  in  a 
good  conductor,  and  the  force  (and  thus  the  voltage)  is  proportional  to  speed.  Modem  Hall 
sensors  use  silicon  for  its  high-resistance  material  characteristic  and  not  for  its  switching  and 
control  capabilities.  A  typical  Hall  sensor  of  the  linear  variety  may  have  a  high  resistance  (about 
1000  W)  across  the  reference  current  leads  and  perhaps  a  little  less  resistance  across  the  sense 
leads. 


Table  1  summarizes  advantages  and  disadvantages  of  the  sense  coil;  table  2  summarizes 
those  for  the  Hall  effect  sensor. 

Table  1.  Advantages  and  Disadvantages  of  the  Sense  Coil 


ADVANTAGES 

DISADVANTAGES 

Extremely  reliable 

No  built-in  electrical  isolation 

No  external  power  requirements 

Comparatively  large  size  and  weight 

Zero-output  for  zero-magnetic  field 
temporal  variation 

Output  amplification  required 

Output  is  a  function  of  magnetic  field 
frequency 

Measures  only  time-varying  magnetic  fields 

Easily  understood  in  the  context  of 

Faraday’s  and  Ohm’s  Laws 

Comparatively  expensive 

3 


Table  2.  Advantages  and  Disadvantages  of  the  Hall  Effect  Sensor 


ADVANTAGES 

DISADVANTAGES 

Built-in  electrical  isolation 

Magnetic  field  strength  and  power  supply 
variations  must  be  accounted  for 

Extremely  small  size  and  weight 

Requires  an  external  power  supply 

Very  reliable 

Offset  for  zero-magnetic  field 

Low  cost 

Less  easily  understood  than  sense  coil 

Measures  both  constant  and  time-varying 
magnetic  fields 

3.  ANALYTICAL  MODELING 


An  analytical  model,  developed  using  the  Biot-Savart  law,  Maxwell's  field  equations,  and 
the  Evershed  model  of  PM  bars,  predicts  the  magnitude  and  frequency  of  the  voltage  waveform 
induced  in  a  copper-wire  sense  coil.  That  model  also  helps  define  the  requirements  of  the  data 
acquisition  and  peripheral  instrumentation  control  electronics  for  the  experiment  prior  to  actual 
high-speed  projectile  tests.  The  expense  of  those  tests  dictate  that  all  issues  regarding  data 
acquisition  and  control  be  resolved  before  the  start  of  testing. 

Integration  of  the  expression  for  the  normal  component  of  a  magnetic  dipole  field  was 
carried  out  over  the  area  enclosed  by  a  sense  coil  loop  for  horizontal  orientation  of  the  moment. 
The  resulting  formula  for  0,  the  magnetic  flux  linked  to  coil  loop,  is  a  function  of  (1)  the 
dimensions  of  the  coil  loop,  (2)  the  distance  from  the  center  of  the  loop  to  the  dipole,  and  (3)  the 
magnitude  and  orientation  of  the  PM  moment.  Rectangular  coordinates  were  used  in  which  z  was 
the  axial  distance  of  the  moving  permanent  magnet  m,  and  r  was  the  radius  of  the  coil  loop  (see 
figure  1). 

A  PM  can  be  represented  by  a  single  magnetic  dipole.  Accordingly,  the  flux  expressions 
were  derived  on  that  basis,  taking  into  account  the  relative  position  of  the  coil  loop  and  the 
position  of  the  moving  PM.  A  procedure  was  outlined  for  computing  0  when  the  PM  dipole  was 
oriented  in  an  arbitrary  manner  or  when  it  was  desired  to  obtain  0  for  a  space-distribution  of 
dipoles. 


The  effects  of  variations  of  strength,  length,  and  diameter  of  the  PM  on  0  were  studied, 
the  results  of  which  were  used  to  determine  the  size  of  the  coil  loop  for  optimum  response. 
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The  magnetic  induction  at  point  P  due  to  the  presence  of  a  small  bar  PM  of  pole  strength 
q  at  the  origin  of  the  coordinate  system  (see  figure  3)  is  given  by 

B  =  fjH  =  (3M*rr  -r2M)r~5 ,  (3) 

where  M  =  ql  is  the  total  magnetic  moment  of  the  PM,  and  /,  its  length,  is  much  less  than  r  (the 
magnitude  of  the  radius  vector  r).  If  the  quantities  on  the  right  side  of  equation  (3)  are  measured 
in  ampere/meter,  B  is  in  teslas. 


Figure  3.  Geometry  Used  in  Developing  the  Magnetic  Flux  Equations 

The  components  of  B  with  respect  to  any  arbitrary  Cartesian  coordinate  system  with 
origin  at  point  D  are  found  by  the  usual  methods  of  reduction  to  be 


Mx{3x 2  -  r2>)  +  3Myxy  +  3Mzxz\ 

r 

(4) 

bAfja:  +  My{3y2  -  r2J  +  3Mzyz\ 

By  5  ) 

(5) 

3 M  zx+  3M  zy  +  M,(3z2  -r2)] 

Bz  =  pHz  =  [- - - -A - 11, 

(6) 
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where,  according  to  figure  3, 

Mx  -  \M\  cos(«), 

My  =  \M\  cos(/?), 

Mz  =  |M|cos(/), 
and 

r 2  —  x1  +y2  +z2 .  (7) 


The  quantity  to  be  studied  is  the  integral  of  Bz  over  a  circular  region  in  an  xy-plane  parallel  to  the 
coordinate  XY-plane  for  special  cases  of  the  angles  a ,  f3,  and  y .  This  integral  predicts  the 
magnetic  field  signatures  from  configurations  of  one  or  more  permanent  magnets. 


3.1  PERMANENT  MAGNET  MODELING 

A  PM  can  be  modeled  as  a  magnetic  dipole  pair,  with  the  strength  of  each  pole  <f>  equal  to 
the  flux  through  the  central  section  of  the  magnet  and  the  pole  spacing  /p  is  almost  equal  to  the 
magnet  length.  The  magnet  is  assumed  to  be  cylindrical  with  square-cut  ends,  axially  magnetized, 
and  made  of  a  low  permeability  material  (such  as  a  ferrite  or  neodymium-boron-iron  (Nd  B  Fe). 

Using  the  Evershed  relationship,  the  demagnetizing  factor*/^,  can  be  calculated  as 


J 'DM  ~  d  ^1  +  £>’ 

where  D  is  the  diameter  of  the  PM  and  X  is  its  length. 

The  magnetic  moment  of  the  test  magnets  were  measured  by  inserting  them  into  the 
working  region  of  a  Helmholtz  Coil  pair,  with  proper  orientation.  The  coil  was  connected  to  an 
integrating  fluxmeter.  The  meter  was  zeroed,  and  the  magnet  was  then  withdrawn  a  considerable 
distance  and  rotated  90°  relative  to  the  axis.  The  magnetic  moment  was  then  found  from  the 
product  of  the  coil  constant  and  fluxmeter  reading,  as  described  in  reference  2. 

In  the  mid- 1800s,  Ampere  and  others  proved  that  steel  bar  magnets  could  be  modeled  as  if 
they  consisted  of  two  poles  of  magnetic  flux,  a  source  (north  or  +)  and  sink  (south  or  -),  and  that 
the  flux  travels  in  or  out  of  the  pole  evenly  in  all  directions  (see  figure  4). 
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Figure  4.  Permanent  Magnet  Model 


In  general,  <fr- 


mmf 
R  ’ 


(9) 


where  mmf  is  the  magnetomotive  force,  R  is  the  reluctance,  and  <j>  is  the  flux.  If 
t  =  \BdA, 


and 


mmf  =  Ji7  •  dl. 


(10) 

(11) 


where  B  is  the  flux  density,  Ax  is  the  cross-sectional  area  of  the  magnet,  and  H  is  the  magnetic 
field  intensity,  it  follows  that 


BA  = 


H-L 


eff 


R 


and 


B  L 


H  AR ’ 


(12) 

(13) 


where  Leff  is  the  effective  distance  between  magnet  poles. 

Set  Lefr  =k„  —  , 

p  AxR 

where  L  is  the  magnet  length  and  kp  is  the  ratio  of  pole  spacing  to  total  length 

±=k,JL. 

H  p  AR 


(14) 


(15) 
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Assuming  the  two  poles  behave  independently,  the  reluctance  in  air  is 


R  =  ±  fA. 

A/ 


where  As  is  a  total  spherical  area.  Because  the  flux  is  uniform  in  all  directions, 
A.  =  4nr2 


(16) 


and 


R 


It '  dr  _  1  r 

J  A-rrr ~  Atth  « 


(jne  j?o/e 


i  n 


ju0  J  4 nr'  4xju0J  4tv~  4;r//0  \aJ  4^//0 


1  l 


(17) 

(18) 


The  upper  limit  of  integration  should  be  set  at  infinity.  If  the  lower  limit  is  set  at  zero,  the 
reluctance  becomes  infinite.  Instead,  assume  a  lower  limit  of  a  sphere  of  area  As  and  radius  rs. 
The  reluctance  of  one  pole  is  then 


R„ 


1 


om  pole  A 

4  XMors 

The  total  reluctance  of  both  poles  is  therefore 

R Total 


1 


2^0  r, 

The  area  of  one  of  the  spheres  is 
A,  =  4  7rrJ 


and 


( Ar\  2 


U  n) 


The  reluctance,  in  terms  of  sphere  area,  A  s ,  is 


R  = 


1 


2 7TjU0  \4n; 

Setting  the  area  of  two  spheres  equal  to  the  actual  magnet  area  As 


R- 


1 


7tA, 


(19) 


(20) 


(21) 

(22) 


(23) 


(24) 
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where  As.  =  2A*  Substituting  this  result  into  the  B/H  formula, 


(1  Y 

B  ,  Mi*1' 

H  =  k ' - T- 


For  a  round  bar  magnet  with  square-cut  ends, 
Y  =  (side  area )  +  2  ■  (end  area). 


and 


A  =tzDL  +  2-  \-Dz  \, 

\4 


Ax  =  —D2. 
x  4 


Substituting  these  results  into  equation  (25), 


-B  4  L 

M0H~  p  nD2 
and 


"l 

r 

^D2  V 

—  n 

kDL  + 

2 

V 

2  Y 

-B 


=  k 

MoH  pKD 


21  v,  2  L)2 


1  + 


D 


(25) 


(26) 

(27) 


(28) 


(29) 


3.2  CALCULATION  OF  MAGNETIC  FLUX  DENSITY 

The  average  flux  density  B  can  be  found  from  the  B/H  curve  of  the  permanent  magnet 
material.  If  the  recoil  permeability  of  the  magnet  material  is  known,  and  if  the  B/H  line  remains  in 
the  linear  region, 


R  —  f? 

Average 


1  + 


Br 


f  DM  * 


(30) 


For  ferrite  and  Nd-B-Fe,  /ur  is  about  1.04  to  1.08,  and  an  average  value  is  fj.r  =  1.06.  The  flux 
(equal  to  the  pole  strength)  is 

$  ~  ^ Average  '  ^ M  >  (^) 
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(32) 


Am=ttD2. 

The  flux  can  be  measured  with  a  search  coil  closely  fitting  the  magnet  and  an  integrating 
fluxmeter.  The  magnetic  moment  m  of  the  magnet  is 

(33) 

where  lp  is  the  spacing  between  equivalent  poles. 


3.3  CALCULATION  OF  INDUCED  VOLTAGE  (BACK  ELECTROMOTIVE  FORCE) 

The  voltage  induced  in  a  pickup  coil  is 

E  =  N^f-,  (34) 

cit 

where  N  is  the  number  of  coil  turns  and  (f>c  is  the  flux  passing  through  the  coil. 


The  flux  linking  a  coil  at  distance  z  from  a  coil  of  radius  r ,  from  a  magnetic  dipole  of 
strength  m  is  calculated  by  assuming  that  at  the  positive  (N)  pole,  the  pole  is  a  flux  source,  with 
flux  <j>  radiating  evenly  in  all  directions.  At  the  opposite  (S)  pole,  a  flux  sink  exists,  and  flux 
flows  inward  from  all  directions  evenly  (the  total  amount  of  which  is  (f> ). 


As  shown  in  reference  2, 

f  \ 

A  1  ^ 

<f>c  =  m - - - r  ■ 

The  voltage  induced  in  the  coil  is  expressed  as 


(35) 


E  =  N 


d <K 

dt 


=  N 


d4>c  f  dz\ 
dz  \dt )’ 


(36) 


where  —  is 
dt 


the  velocity  of  the  magnet. 


It  follows  that 


r 


dK 

dz 


-  3  m 


V 


A 


(r2+z>) 


(37) 
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and 


E(N,r,  z,  v,  L,D,B„m,)  =  N^-, 


3  Nm 


(r2  +  z2 


)3  (r‘+z’)\ 


dz 

Tt 


(38) 


where  N  is  the  number  of  turns  in  a  pickup  coil,  and  r  is  its  radius;  z  is  the  axial  position  of  the 
PM,  L  is  its  length,  D  is  its  diameter,  m  is  its  moment,  and  jj.r  is  its  permeability;  and  Br  is  the 
remanent  flux  density. 

Figure  5  shows  an  analytical  model  that  calculates  the  EMF  of  a  PM  falling  under  the 
influence  of  gravity  using  the  previously  derived  equations  and  data.  A  plot  of  generated  voltage 
for  the  parameters  given  (table  3)  is  shown  in  figure  6  and  corresponds  well  with  the  experimental 
data  shown  in  figure  7.  Figure  8  superimposes  the  experimentally  measured  and  calculated 
voltages  in  single  plot.  Figure  9  re-plots  the  analytical  data  with  respect  to  position. 
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Figure  5.  Analytical  Model  of  a  Falling  PM  (Voltage  Equation) 
Table  3.  Parameters  To  Determine  Generated  Voltage 


N  :=  150 

Number  of  coil  turns 

pi  :=  3.14159265 

r :=  .23495 

Radius  of  coil 

hi  :=-i.o 

Position  1 

h2  :=  1.0 

Position  2 

g  :  =  9.807 

Gravitation 

z  :=  hi, hi  +  0.01..  h2 

d  :=(■!)■£* 

\4/  100 

Diameter  of 
permanent  magnet 

Bavg  :  -  '  1.1564 

Permanent  magnet 
flux  density 

d  =  0.006 

2.54 

1  •=  (0.377) - 

100 

Length  of 
permanent  magnet 

1  =  0.01 

A  :=£?•<? 

Permanent  magnet 

A  =  3.167- 10-5 

4 

cross-sectional  area 

v  :=  A-i 

Permanent  magnet 
volume 

V  =  3.033*10  7 

M  • =  Bavg-  V 

Magnet  moment 

M  =  -3.507*  10  7 
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Figure  6.  Generated  Voltage  as  a  Function  of  Time 


Figure  7.  Experimentally  Measured  Voltage  as  a  Function  of  Time 


4.  FEA  MODELING 


A  2-D,  axisymmetric  FEA  model  was  created  and  analyzed  to  calculate  the  magnetic  flux 
distribution  and  to  verify  the  analytical  model  and  experimental  measurements  (references  3  and 

4). 

Figures  10  and  1 1  show  the  axisymmetric  2-D  FEA  model.  The  PM  permeability  and 
coercive  force  were  obtained  from  the  magnet  distributor.  The  volume  surrounding  the  PM 
measures  40  inches  in  the  axial  direction  and  30  inches  in  the  radial  direction. 


Figure  1 0.  2-D  Axisymmetric  FEA  Model 


Figure  11.  Expanded  View  of  the  2-D  Axisymmetric  FEA  Model 
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Figure  12  shows  the  magnetic  flux  lines  in  the  surrounding  volume  that  is  shown  in  figure  1 1 . 
Figures  13  and  14  illustrate  the  magnetic  flux  density  in  the  axial  and  radial  directions  as  a 
function  of  the  radial  distance  from  the  PM.  These  FEA  magnetic  flux  density  calculations  are 
compared  with  the  analytical  model  and  show  close  correlations. 


Figure  13.  Magnetic  Flux  Density  in  the  Axial  Direction 
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Figure  14,  Magnetic  Flux  Density  in  the  Radial  Direction 


4.1  PERMANENT  MAGNET  MATERIAL  SELECTION  USING  AN  FEA  MODEL 

An  FEA  model  was  built  and  analyzed  to  find  the  effect  of  PM  material  in  the  projectile. 
Figures  15  and  16  show  the  different  results  from  two  different  PM  strengths.  The  results  in 
figure  15  from  the  40  mega-gauss  oersted  (MGOe)  with  979  IcA/m  Nd  B  Fe  PM  and  the  magnetic 
flux  density  are  plotted  as  a  function  of  the  PM  position.  Figure  1 6  shows  the  results  from  the  22 
MGOe,  764  kA/m  PM  material,  revealing  that  the  stronger  PM  produces  a  15-percent  higher 
magnetic  field  in  the  air  around  it.  From  those  plots  the  maximum  magnetic  flux  density  was 
calculated. 
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Figure  15.  Magnetic  Flux  Density  as  a  Function  of  Magnet  Position  (40  MGOe,  979  kA/m) 
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Figure  16.  Magnetic  Flux  Density  as  a  Function  of  Magnet  Position  (22  MGOe,  764  kA/m) 


4.2  PROJECTILE  STEEL  MATERIAL  SELECTION  USING  AN  FEA  MODEL 

The  FEA  model  was  also  used  to  find  the  appropriate  material  for  the  projectile. 
Permeability  of  the  casing  surrounding  the  PM  was  varied  in  appropriate  intervals  from  1.0  to 
1 50.  Higher  permeability  means  that  the  magnetic  flux  encounters  less  reluctance;  therefore, 
more  flux  flows  in  the  material.  In  electric  motor  design,  higher  permeability  materials  are  always 
preferable  for  flux  concentration.  In  this  project,  however,  a  low  permeability  material  is  desired 
so  that  the  flux  is  allowed  to  flow  outside  the  projectile. 

Figures  1 7  and  1 8  show  the  magnetic  flux  lines  of  the  same  PM  strength  but  two  different 
casings  with  permeabilities  of  1 .0  and  10.0  respectively.  From  these  results,  it  was  determined 
that  permeability  should  be  less  than  5  to  assure  the  broadest  possible  measurement  range.  Levels 
as  high  as  1 0  are  still  practical  depending  on  the  surrounding  magnetic  noise  level. 
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Figure  1 7.  Flux  Distribution  Plot  fa- 1. 0) 


Figure  18.  Flux  Distribution  Plot  fa=  10. 0) 


5.  EXPERIMENTAL  MEASUREMENTS 


5.1  SENSE  COIL  MEASUREMENTS 

Experiments  with  the  sense  coil  were  conducted  by  first  dropping  a  PM  through  the  coil 
center  and  then  shooting  a  PM  projectile  horizontally  through  the  center  of  the  coil  using  an  air 
compressor.  In  both  cases,  the  PM  moved  through  a  smooth  glass  tube.  The  analytical  and  FEA 
models  developed  proved  to  be  an  accurate  means  of  predicting  the  induced  voltage.  From  this 
calculated  induced  voltage,  the  position  and  speed  of  the  projectile  were  determined. 

Figure  19  shows  a  block  diagram  of  the  equipment  used;  figure  20  shows  the  test  setup 
picture  for  the  permanent  magnet  drop  test.  Several  PM  drop  tests  were  performed  and  induced 
coil  voltage  waveforms  recorded.  Combinations  of  1,  3,  and  7  PMs  attached  in  series  were  used 
for  this  experiment.  Each  PM  was  made  of  Nd  B  Fe  and  was  cylindrical  with  dimensions  of  a 
0.249-inch  diameter  and  a  0.377-inch  in  length.  Magnet  poles  were  located  at  the  top  and  bottom 
of  the  cylinders. 

The  results  of  a  drop  test  using  one  PM  is  shown  in  figure  7.  Measurements  from  tests 
with  three  and  seven  PM  combinations  are  shown  in  figures  21  and  22.  Figure  23  shows  the 
shout-test  setup  picture.  The  simulated  and  measured  voltages  were  compared  and  showed  about 
a  3 -percent  error  between  them. 


Figure  19.  Permanent  Magnet  Drop-Test  Setup  Block  Diagram 
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Figure  20,  Perman  eni  Magnet  Drop-  Test  Setup  Picture 
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Figure  22.  Drop- Test  Measurements  with  Seven  Permanent  Magnets 


Figure  23.  Shout-Test  Measurements  Test  Setup 


5.2  REPEATABILITY  AND  ACCURACY  OF  MEASUREMENTS 


Experimental  test  data  it  must  be  accurate  and  repeatable  to  be  useful.  A  test  was 
performed  whereby  the  same  set  of  data  was  recorded  twice  and  the  difference  between  the  two 
was  computed.  The  calculated  differences  were  close  to  zero.  The  test  was  executed  at  several 
different  speeds,  and  each  time  the  results  were  nearly  the  same. 


5.3  HALL  EFFECT  SENSOR  MEASUREMENTS 

Static  flux  density  measurements  were  performed  using  a  commercially  available  Hall 
effect  sensor.  Those  measurements  confirmed  the  manufacturer’s  claim  that  voltage  varied 
linearly  with  flux  density  (figure  24).  Because  the  sense  coil  experiments  verified  the  results  of 
the  analytical  and  numerical  models,  no  further  experiments  with  moving  PMs  were  performed. 


Figure  24.  Voltage  Varies  Linearly  with  Flux  Density 
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6.  CONCLUSION 


Analytical,  numerical,  and  experimental  techniques  were  combined  to  solve  the  problem  of 
accurately  determining  the  position  and  velocity  of  a  high-speed  projectile.  First,  an  analytical 
model  was  created  to  predict  the  behavior  of  falling  (low-speed)  permanent  magnets.  This  model 
and  a  more  sophisticated  FEA  model  were  used  to  select  appropriate  materials  and  sensing 
equipment  for  running  experiments.  These  experiments  verified  the  analytical  and  numerical 
models,  proving  that  the  analytical  and  numerical  models  were  also  useful  for  sensing  the  behavior 
of  high-speed  projectiles. 


REFERENCES 


1.  P.  Emerard  and  J.  Gilbert,  “Introduction  Theory  and  Fundamental  of  Hall  Effect  Sensor  and 
Transducer,”  IPCIM ,  September  14,  1995. 

2.  J.  Stupak,  “A  Method  of  Calibrating  Helmholtz  Coils  for  the  Measurement  of  Permanent 
Magnets,”  24th  Annual  Symposium,  IMCSD,  June  1995. 

3.  D.  Lowther  and  P.  Silvester,  Computer-Aided  Design  in  Magnetics ,  Springer-Verlag, 
Berlin,  New  York,  1986. 

4.  M.  Yaksh,  “Electromagnetic  Analysis  of  a  Machine,”  ABB  Impell  Corporation,  (in 
preparation). 


23/24 

Reverse  Blank 


NUWC  NPT  Technical  Memorandum  962063 


External 

Internal 

Codes: 


DISTRIBUTION  LIST 


01 

0212 

51122  (routing  copy) 

5131  (Conforti) 

5141  (2)  • 

10 

102 

1 04  (Ke  shura) 

20 

22 

34 

38 

382 

383 
60 
80 
81 
814 
82 
821 
822 
8222 

8224 

8225 
823 

823 1  (Barnett,  Cho  (5),  Egan,  Raposa,  Thivierge,  Atwater) 

8232 

8233 

8234 

8291  (Benson,  Cancilliere,  Hillenbrand) 

8292 
83 

832 

8322  (Curtis  (1),  Kirschner  (1)) 

833 


Total:  48 


